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We report the discovery of superconductivity below 1.65(6) K in Sr-intercalated graphite SrCe, 
by susceptibihty and specific heat (Cp) measurements. In comparison with CaCs, we found that the 
anisotropy of the upper critical fields for SrCe is much reduced. The Cp anomaly at Tc is smaller 
than the BCS prediction indicating an anisotropic superconducting gap for SrCe similar to CaCe- 
The significantly lower of SrCg as compared to CaCe can be understood in terms of "negative" 
pressure effects, which decreases the electron-phonon coupling for hath in-plane intercalant and the 
out-of-plane C phonon modes. We observed no superconductivity for BaCg down to 0.3 K. 

PACS numbers: 74.70. Ad, 74.25.Bt, 74.62.Fj, 74.25.Kc 



The discovery of superconductivity in YbCe and CaCe 
P, Q initiated intensive theoretical and experimental 
investigations [3, 4, i, H, 0, IE S E Ell, [ll on the 
alkahne-earth graphite intercalation compounds (GICs). 
The results of the experimental studies, specially the 
observation of a Ca isotope effect strongly favor 

electron-phonon(e-ph) coupling rather than exotic elec- 
tronic mechanisms j^. Recent ab-initio electronic struc- 
ture calculations prove that, in contrast to MgB2, e-ph 
coupling involving electronic intcrlayer (IL) states be- 
comes relevant [3] and is sufficiently strong to generate 
the relatively high Tc's 4,, 9]. These findings stimu- 
lated discussions on the variability of the e-ph coupling 
strength in different branches of the electronic and vi- 
brational states for honeycomb layered compounds like 
MgB2, CaSiJlll, (Ca,Sr)AlSi [li| and also the hypothet- 
ical Li2B2 [iq . 

Many issues still remain open, especially about the na- 
ture of the relevant phonons. In order to modify the 
phonon spectrum, Ca isotope substitution [13] and hy- 
drostatic pressure experiments [1, [3, have been per- 
formed. The isotope experiments reported a surprisingly 
high isotope exponent for Ca, Q:(Ca) « 0.5, close to the 
BCS limit, suggestin g a dominant role of the Ca phonons 
in the e-ph coupling |12l |. However, ab-initio calculations 
predicted similar isotope exponents ~ 0.25 for Ca and 
C, pointing to comparable contributions to the e-ph cou- 
pling from the Ca and C phonons [1] . The positive pres- 
sure dependence of Tc found for CaCg has been discussed 
in terms of phonon softening for the in-plane Ca vibra- 
tions 0, Although the experimentally determined 
Tc's grow almost linearly with pressure, the ab-initio cal- 
culations predicted a non-linear increase with a reduced 
magnitude [§]. To resolve these discrepancies, possible 
anharmonic effects of the ultrasoft intercalant phonon 
modes, or a continuous superconducting gap distribu- 
tion due to anisotropic e-ph coupling have been suggested 



17|, [18|, asking for further investigations. 
Another way to modify the relevant phonon modes is 



to vary the intercalant species and replace Ca with other 
alkaline-earths such as Sr or Ba. Mazin pointed out that 
for CaCe and YbCe the square root of the mass ratio 
of the intercalants is only 15% larger than the ratio of 
their Tc's Thus, according to this "isotope" effect 
argument other alkaline-earths GICs may as well be su- 
perconducting. In fact, subsequent ab-initio calculations 
predicted superconductivity for SrCg at 3.1 K and BaCg 
at 0.2 k[16|. In this Letter, we report the discovery of 
superconductivity in SrCg at Tc — 1.65(6) K by suscep- 
tibility and specific heat measurements and the absence 
of superconductivity in BaCg down to ^ 0.3 K. The su- 
perconducting properties of SrCg as well as the ab-initio 
calculations clearly demonstrate that SrCg can serve to 
bridge the two seemingly different classes of the super- 
conducting GICs: the low- Tc alkali GICs and the newly- 
discovered " high- Tc" systems, CaCe or YbCe. Further- 
more, the comparison of SrCg with CaCg provides a bet- 
ter insight into the unconventional nature of supercon- 
ductivity in alkaline-earth GICs. 

Samples of SrCg and BaCe were synthesized from 
pieces of highly oriented pyrolytic graphite (Advanced 
Ceramics, size ~ 3x1x1 mm"^) and Sr (99.95%) or Ba 
(99.95%) metal by a vapor phase reaction performed for 
more than a month at 470°C and 500°C, respectively. X- 
ray diffraction patterns show no reflections due to pristine 
graphite or other higher stage intercalated phases, indi- 
cating good sample homogeneity. The graphite layer dis- 
tance increases from CaCg (4.50 A) to SrCg (4.95 A) and 
to BaCg (5.25 A) as expected from the size of the inter- 
calant atoms. The stacking sequence of SrCg and BaCfiis 
found to be of the a/S-type (space group PGs/mmc) [l9| 
which differs from the a/37- type stacking in CaCg- Super- 
conductivity for SrCe was first observed by ac magnetic 
susceptibility (119 Hz) measurements to 0.3 K, and sub- 
sequently confirmed by specific heat measurements using 
a PPMS '^Hc calorimeter (Quantum Design). The pres- 
sure dependence of Tc was measured up to ~ 1 GPa as de- 
scribed in detail elsewhere Ml . In order to understand the 
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superconducting properties of SrCe and compare with 
those of CaCe, we also performed ab-initio calculations 
of the electron-phonon properties for SrCg using the ex- 
perimental a/? stacking 20l. |21|. 



Figure 1 shows a sharp superconducting transition in 
the magnetic susceptibility x(r) at ~ 1.65 (6) K in 
SrCg. The transition width A 0.06 K was deter- 
mined as the temperature difference between 10% and 
90% of the diamagnetic shielding. In a magnetic field, 
Tc shifts to lower temperatures and the superconduct- 
ing transition broadens. In BaCe, we cannot find any 
signature of superconductivity down to 0.3 K. 

The upper critical fields, iJc2^' ^'^^ ^ parallel and per- 
pedicular to the c axis follow the Werthamer-Helfand- 
Hohenberg (WHH) prediction rather well ?233(Fig.[l]^b)). 
A slight deviation from the WHH curve is found at low 
temperatures, which has been observed in other GICs 
such as CaCe and KCg [2^. Within the scope of the 
WHH approximation, we obtained ffj.'2(0) ~ 138 Oe and 
H^{0) w 276 Oe as well as the corresponding coherence 
lengths, Ca6(0) « 1510 A and Cc(0) w 700 A. Similar to 
CaCe, Cc(0) is larger than the graphite layer distance, 
indicating 3-dimensional (3D) superconductivity. 

Bulk superconductivity in SrCg is confirmed by specific 
heat (Cp) measurements. The characteristic Cp anomaly 
at Tc — 1.65 K is clearly observed at H = 0, and com- 
pletely disappears with H = 500 Oc {H > i?c2(0)). 
There is no offset of Cp/T at as T ^ K, in- 
dicating a complete superconducting phase. Similar to 
CaCg, the normal state Cp deviates slightly from a 
dependence (the inset of Fig. [2]), due to low- lying Ein- 
stein phonon modes. The normal state C„ can well be 




FIG. 1: (Color online) (a) Temperature dependence of the ac 
susceptibility of SrCe for various fields (as indicated) perpen- 
dicular to the c-axis. Tc is determined from the intersection 
between the extrapolated lines of the steepest slope of xiT) 
and of the normal state x{T) as shown by the solid (black) 
lines, (b) The Hc2{T) for H \\ ah plane and H \\ c axis. The 
WHH curve for both H directions are shown with (red) lines. 



described by Cp{T) = + Cph{T) where 7jv is the 

Sommerfeld constant, and Cph{T) = (iT^ + ST^ is the 
lattice contribution. The best fit to the H — 500 Oe data 
yields -/n = 5.92(1) mJ/mol K^, f3 = 0.191(1) mJ/mol 
K'^, and S = 0.801(4) ^J/mol K^^. The estimated Dc- 
bye temperature O_d(0) = 414(1) K is lower than that of 
CaCe (0_d(O) = 598 K) as expected from the atomic mass 
difference. From a comparison of with the calculated 
density of states at the Fermi level Ep, N{0), we estimate 
the e-ph coupling strength A using the relation = 
(27r2fc|/3)7V(0)(l-t-A). With N{0) = 1.63 states/eV-cell, 
we arrive at A = 0.54(1), somewhat lower than in CaCe, 
but still in the intermediate coupling regime. 

The difference ACp between the normal and the su- 
perconducting state is shown in Fig. [21 At low temper- 
atures ACp(r) exceeds the BCS prediction while it is 
slightly lower than the BCS value near T^. Using the 
'a-model' which assumes an isotropic s-wave BCS gap 
A(T) scaled by the factor a = A{0) /ksTc, we were able 
to fit the detailed temperature dependence of ACp{T)/T 
by adjusting the gap ratio to a = 1.67. Accordingly, the 
corresponding 2A(0)/A:Brc = 3.34 is somewhat reduced 
below the weak coupling BCS limit of 3.52. 

First, we discuss the anisotropy of the superconducting 
properties of SrCe. The anisotropy of the upper critical 
fields, Th = H^2/h12 amounts to w 2 at T Tc/2, close 
to that of YbCg [IJ , but significantly smaller than found 
in CaCg {Th ~ 4) Assuming an isotropic supercon- 
ducting gap, the anisotropy Th reflects the anisotropy 
of the Fermi velocities. Our ab-initio calculations for the 
Fermi surface (FS) of CaCe and SrCg clearly reveal an el- 
liptical sheet, associated mainly to interlayer states, and 
a tubular structure, of tt* character (c/. Ref. [lB|). In 
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FIG. 2: (Color online) Temperature dependence of ACp/T 
= Cp{H^ 0)/T - Cp{H = 500 Oe)/r. The (blue) dashed 
and (red) solid lines represent the BCS curve and the best fit 
according to the a-model (see text), respectively. The inset 
shows the temperature dependence of Cp at H— and 500 
Oe. The solid (black) line through the data points for H = 
500 Oe is a fit to a polynomial as described in the text. 
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SrCg the FS for the IL bands has a more pronounced 
2D character and is open, which results from the 10% 
increase of the c axis lattice parameter as compared to 
CaCg- The anisotropy of the average Fermi velocities 

1/2 



Af = (< Vp^^ >/< is very close in the two com- 



pounds: Fyp « 1.9 in CaCe and f» 1.7 in SrCe. If we, how- 
ever, consider only the IL sheets where the superconduct- 
ing gap is larger [l^, we find a much larger difference: 



1.1 for CaCe and = 2.1 for SrCg. Therefore, 



Fh of SrCg is expected to be larger or at least similar 
to that of CaCg. The significantly enlarged anisotropy 
of Hc2 in CaCe with respect to SrCg therefore cannot 
be understood in terms of the anisotropy of the Fermi 
velocities but must be attributed to an anisotropy of the 
superconducting gap as well. 

This conclusion is supported by the reduced Cp jump 
observed in SrCg. Assuming an isotropic gap, the nor- 
malized jump, ACp{Tc)/jNTc, grows with the e-ph cou- 
pling strength over the BCS weak coupling limit of 1.432. 
With A w 0.54 estimated from Cp, ACp{Tc)/jnTc is 
expected to be enhanced over the BCS value. The ex- 
perimental ACp{Tc)/jNTc = 1.426, however, is smaller 
than the BCS value. As a characteristic feature of the 
anisotropic superconducting gap, the entropy "lost" near 
Tc is transferred to lower temperatures (Fig. [2]) [l^ . 
These findings indicate that the superconducting gap in 
CaCe as well as in SrCe has a marked anisotropy, which 
is also supported by recent calculations [3] ) showing that 
CaCe indeed exhibits a strongly k-dependent supercon- 
ducting gap due to anisotropic e-ph interaction. Consid- 
ering that the deviation of ACp(T) from the predicted 
curve for the isotropic gap model is less pronounced 




FIG. 3: Tc as a function of the graphite layer distance, d for 
the alkali-GICs, XCs {X — K, Rb and Cs) and the alkaline 
earth-GICs XCq {X = Ca, Yb^r, and Ba). For CaCe, Tc 
at high pressure (P = 8 GPa)[l0| is also plotted (the half- 
shaded square) and the graphite layer distance for the com- 
pressed CaCg is estimated from the theoretically calculated 
bulk modulus 0]. The upper limit of Tc for BaCs is indicated 
by the arrow. The inset shows Tc vs. pressure for SrCe. 



for SrCg than for CaCe fv^, the superconducting gap 
anisotropy is weaker for SrCg. This is also consistent 
with the reduced anisotropy in Hc2- Therefore, replacing 
Ca with Sr decreases not only the strength of the e-ph 
coupling, but also its anisotropy. 

We can now discuss the reduced Tc of SrCg, which is 
about an order of magnitude smaller than of CaCg. 
When Ca is replaced with Sr, Tc decreases much more 
than it would be expected in the view of the alkaline- 
earth "isotopic" substitution i. With ^A/sr/Mca 0^ 
1.48, one would expect a Tc of 7.8 K for SrCe, hence 
much larger than the experimental value. We thus con- 
clude that the mass of the intercalant is not the main 
factor which determines Tc. Instead, we found that the 
Tc's of the superconducting GICs strongly depend on the 
graphite interlayer distance, d. Figure [3] illustrates that 
Tc decreases rapidly, almost exponentially with increas- 
ing d for both alkali- and alkaline earth-GICs [25, 261 . 
The increase of Tc with pressure found for CaCe [a Isl) llO| 
and SrCg [dT^/dP ~ 0.35 K/GPa, the inset of Fig. [S]) 
clearly manifests a similar trend. Obviously, the main 
factor that governs the significant decrease of Tc in SrCe 
as well as the absence of superconductivity in BaCe down 
to 0.3 K, is the increased distance between the graphite 
layers. 

Figure 3] clearly demonstrates why. In the calcu- 
lated phonon density of states (PhDOS) (Fig. [Kc)) 
and the corresponding Eliashberg functions a^T(w) 
(Fig. m^d)) for CaCe and SrCe, we observe three groups 
of phonons: intercalant-related vibrations {Ixy and F) 
at a; < 20 meV, C out-of-plane vibrations {Cz) around 
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FIG. 4: (Color online) (a) Tc vs. the position of the bot- 
tom of the interlayer bands {Eb) with respect to Ef (b) Tc 
vs. the band width of the interlayer bands along the c axis 
{Wfi^). The grey line is meant as guide to the eye. (c) 
Phonon density of states and (d) Eliashberg function o?F{Q,) 
and frequency-dependent electron-phonon coupling \{uj) = 
2 a'^F{Q.)/Q.dVl for CaCe (orange) and SrCe (black). 
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50 meV, and C bond-stretching vibrations (Cxy) at 
Lu > 150 meV. The quaUtative shape of a'^F{uj) is the 
same in the two compounds, indicating a similar spec- 
tral distribution of the e-ph interaction. But the to- 
tal A decreases from 0.83 in CaCe to 0.56 in SrCe, 
while the logarithmic-averaged phonon frequency (ujin) 
remains unchanged at ^ 305 K. The results are in very 
good agreement with the previous calculations for SrCg 
with the a/37 stacking Within computational accu- 
racy, the difference between the a/3 and a/?7 stackings 
seems negligible in contrast to a previous conjecture [J]. 
Using the Allen-Dynes formula, with /i* — 0.145, we ob- 
tain Tc = 11.4 K for CaCe and = 3.1 K for SrCg, in 
reasonable agreement with the experimental results and 
previous calculations [l^. 

The reduction of Tc in SrCe is due to the simulta- 
neous decrease of the I^y and Cz contribution to the 
e-ph coupling. The reduced coupling for the low-lying 
Ixy vibrations has a negative effect on Tc, but it is also 
very effective in increasing ojin, thus, the effect is partly 
compensated. On the other hand, the large reduction of 
coupling associated to Cz vibrations, which happens at 
energy scales comparable to tvin, reduces A leaving ujin 
unchanged, and thus has a large effect on the Tc- The re- 
duction of e-ph coupling for Cz vibrations in the interme- 
diate energy range is essential to explain the significantly 
lower Tc for SrCe than CaCg. 

The e-ph interaction for the Cz modes, and to a lesser 
extent for the Ixy ones, is associated to the interband 
couphng between the IL and tt* states (i^. Going from 
CaCg to SrCg, this interband coupling is essentially re- 
duced by the increase of the c axis lattice constant, which 
decreases the real-space overlap between IL and tt* wave- 
functions. The characteristic parameter monitoring the 
IL-TT* overlap is the bandwidth of the IL band along 
the c axis {Wfj^). Figure IHb) demonstrates that Tc for 
both alkali- and alkaline earth-GICs is directly corre- 
lated to Wf^, clearly explaining the Tc dependence on 
the graphite layer distance. This observation is in strong 
contrast to the dependence of Tc on the degree of filling 
for the IL band (Fig. [Ja)). Even though this band has 
to be occupied for superconductivity to occur Q , there is 
no clear correlation of Tc with the number of IL electrons. 

Our findings indicate that C phonon modes play a non- 
trivial role in the superconductivity of the GICs, in con- 
trast with the conjecture from the Ca isotope experi- 
ments Il2| . Since C would give a sizable contribution 



to the isotope effect, the total isotope exponent would 
exceed the BCS limit ax^ = 0.5, which requires further 
studies on the effects of C isotope substitution. 

In conclusion, we reported that SrCe becomes super- 
conducting at Tc — 1.65(6) K, but BaCg stays normal 
conducting down to 0.3 K. The reduced Cp jump at Tc 
in SrCe and the reduced anisotropy Hc2 in comparison 
with the corresponding data for CaCg strongly s upp orts 
the idea of an anisotropic superconducting gap 17, 3|. 



We give clear evidence that Tc of the GICs essentially 
depends on the graphite layer distance due to sensitive 
change of e-ph coupling for both in-plane intercalant and 
the out-of-plane C phonon modes. Our results suggest 
that a possible route to increase Tc is to replace Ca by 
smaller atoms, such as Mg. Several attempts to prepare 
pure MgCg have failed so far, but partial substitution of 
Mg or Li, as long as the filling of the IL bands is kept, 
could be possible to reduce d, stabilize the structure, and 
as a result, increase Tc- 
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